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ABSTRACT 

Gas undergoing Bondi accretion onto a supermassive black hole (SMBH) becomes hotter toward 
smaller radii. We searched for this signature with a Chandra observation of the hot gas in NGC 
3115, which optical observations show has a very massive SMBH. Our analysis suggests that we are 
resolving, for the first time, the accretion flow within the Bondi radius of an SMBH. We show that 
the temperature is rising toward the galaxy center as expected in all accretion models in which the 
black hole is gravitationally capturing the ambient gas. There is no hard central point source that 
could cause such an apparent rise in temperature. The data support that the Bondi radius is at about 
4"-5" (188-235 pc), suggesting an SMBH of 2 x 10 9 M© that is consistent with the upper end of the 
optical results. The density profile within the Bondi radius has a power-law index of 1.03^0 21 which is 
consistent with gas in transition from the ambient medium and the accretion flow. The accretion rate 
at the Bondi radius is determined to be Mb — 2.2 x 10~ 2 M Q yr . Thus, the accretion luminosity 
with 10% radiative efficiency at the Bondi radius (10 44 erg s _1 ) is about six orders of magnitude 
higher than the upper limit of the X-ray luminosity of the nucleus. 

Subject headings: accretion, accretion disks - black hole physics - galaxies: elliptical and lenticular, 
cD - galaxies: individual (NGC 3115) - galaxies: nuclei - X-rays: galaxies 



1. INTRODUCTION 

Understanding accretion onto black holes remains one 
of the most active areas of research in astrophysics to- 
day, both for probing black hole properties and because 
of their impact on larger-scale problems in galaxy and 
structure formation. In many nearby galaxies, the model 
of a central black hole gravitationa lly capturing ambient 
gas (Bondi accretion: iBondil 119521 ) is a vital part of our 
understanding of how black holes are accreting. 

The key to understanding the dynamics of gas in the 
systems of interest lies in correctly modeling the behav- 
ior of the accreting gas once it falls within the grav- 
itational influence of the black hole, the Bondi radius 
Rb = 2GMbh/Cs, where Mbh is the mass of the black 
hole and c s is the sound speed of the gas in the vicinity 
of Rb- In the absence of angular momentum, the black 
hole is predicted to gravitationally capture the ambient 
interstellar medium (ISM) surrounding it at the Bondi 

rate Mb — Air\R 2 B pc S: where A = 0.25 for an adiabatic 
process and p is the density of the gas at Rb- For ob- 
served central p values of many nearby elliptical galaxies 
or Sgr A* , observed luminosities are orders of magnitude 
smaller than what would be predicted based on the Bondi 
accretion rat e and assuming a sta ndard 10% radiative ef- 
ficiency (e.g.. iFabian fc Rccs 1995]). These observational 
results imply that the radiative efficiency of the infalling 
gas is exceedingly low and/or that much less material is 
accreted by the black hole than the Bondi formula im- 
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plies. 

The very weak radiation from most nearby massive 
black holes has prompted a significant theoretical 
effort aimed at explaining the very low radiative ef- 
ficiencies and/or accretion rates. The central idea is 
that accretion proceeds via a hot radiativcly ineffi- 
cient accretion flow (RIAF) rather than a canonical 
thin accretion disk. Early RIAF models include 
the advection - dominated accretion flow (ADAF ; 
HchimanJ fT977t IRees et al.l [TMl INaravan fc Yl Hflfll . 
the convectio n- domin at ed accretion flow ( CDAF ; 
INaravan et all l2000t iQuataert fc Gruzinovl l2000f l. 
and the advection-dominated inflow- outflow solution 
(ADIOS; Bla ndford fc Beeel man 1 fl999h . 

Determining which (if any) of these scenarios describes 
low-ix black hole systems is of fundamental impor- 
tance to our understanding of accretion physics and 
black hole demography. Observational work has focused 
on using spatially unresolved spectral information to 
constrain theoretical models. While such studies have 
been successful in ruling out c lassical ADAF models 
in some instances (e .g., Sgr A*; iDi Matteo et al.l 120011 : 
iBaganoff et"al~ll2003f) . the main limitation has been the 
inability of even Chandra to resolve the accretion flow 
inside Rb and directly determine the density profiles of 
the accretion flow, as it is the shape of the density profile 
that most strongly distinguishes the theoretical models. 
While ADAFs, CDAFs, and ADIOS models all predict 
that the compression of gas within the accretion flow 
within the Bondi radius leads to a T(R) oc relation 
for the gas, the models predicted significantly different 
radial density profiles (Section [5]). 

In this Letter, we present for the first time spatially re- 
solved temperature and density profiles within the Bondi 
radius of a supermassive black hole (SMBH). Dynamical 
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Fig. 1. — Smoothed 0.3-6.0 keV C handra image of NGC3115 
with D25 (dc Vaucoulcurs ct al. 1991) optical contour. The inset 
is the 2.0-6.0 keV image of the central 25" X 25" region. There is 
not a point source at the center, but a plateau in the diffuse X-ray 
surface brightness. 

modeling of the stellar kinematics of the bulge of NGC 
3115 has confirmed that it harbors a 1-2 x lCr Mp, black 
hole (jKormendv et all 119961: lEmsellem et al.l 119990 . At 
a distance of 9.7 Mpc (jTonrv et al.l I2001Q . this makes 
NGC 3115 the nearest > 1O 9 M black hole. Despite 
its large mass, the black hole is very radio quiet, with 
a core 5 GHz luminosity of le ss than 4 x 10 25 erg s _1 
Hz -1 ( Franccs chini et allll998l) . several orders of magni- 
tude lower than other black holes of similar mass. The 
black hole is also undetected in X-rays, with a Chandra 
upper limit of less than 10 38 erg s _1 (Diehl & Statler 
2008; K.-W. Wong et al. 2011, in preparation). The es- 
timated Bondi accretion rate based on the previ ous ob- 
served density of w 0.1 cm" 3 (|David et alj 120061) in the 
central part of the galaxy is ~ 0.1 M@ yr" 1 , implying an 
expected luminosity ~ 10 45 erg s -1 if accretion occurred 
with 10% efficiency The lack of radio power argues that 
it is unlikely that jets are interacting strongly with the 
gas inside Rb- 

2. X-RAY OBSERVATIONS AND SPECTRAL ANALYSIS 

NGC 3115 (Figure [TJ was observed with Chandra on 
2001 June 14, 2010 January 27, and 2010 January 29 for 
37, 41, and 77 ks, respectively. All the data were re- 
processed using the Chandra Interactive Analysis of Ob- 
servations (CIAO) software version 4.3 and the Chandra 
Calibration Database (CALDB) version 4.4.1. We used 
a local background that is extracted from the 70"-90" 
annulus centered at the X-ray peak. The background 
region is far enough from the center so that the source 
removed surface brightness of the X-ray emission is basi- 
cally flat beyond 70". The background only contributes 
less than about 1% of the emission (0.5-6.0 keV) within 
the Bondi radius (< 5") and less than about 7% for the 
ambient region beyond the Bondi flow (4"-10"). Errors 
are given at la level unless otherwise specified. 

We extracted spectra in annuli centered at the central 
peak of the extended X-ray emission. This peak was de- 
tected as a source by CIAO wavdetect. However, further- 



analysis shows that a modeled circular Gaussian width of 
this peak emission is twice extended as that of a nearby 
bright point source and that the central source is static 
(varied less than 2a) within and between the observa- 
tions. Thus, it is dominated by a concentration of hot 
gas and unresolved binaries rather than a single compact 
accreting object. We included this peak emission but ex- 
cluded all the other sources detected by wavdetect during 
the analysis. 

The unresolved X-ray emission is mainly contributed 
by unresolved low- mass X-ray binaries (LMXBs), 
cataclysmic variables and coronally active binaries 
(CV/ABs), and the diffuse hot gas component that we 
are interested in. The very soft emission from the gas and 
the very hard emission from the unresolved LMXBs are 
easy to separate spectrally. No appreciable emission is 
expected from the gas above 2 keV, so we can use the 2-6 
keV emission as a lever arm for estimating the amount of 
unresolved LMXB flux in the 0.5-2.0 keV band and sub- 
tract it accordingly if we assume an appro priate spectral 
model for the LMXBs. [I rwin et ail (|2003f) have shown 
that the summed spectra of resolved \ow-Lx (< 10 37 erg 
s" 1 ) LMXBs in the bulge of M31 are very similar to more 
luminous LMXBs. Hence, it is valid to assume that the 
unresolved LMXB emission can be spectrally modeled as 
the brighter resolved sources. We found that the fainter, 
softer source like those that make up the Galactic Ridge 
emission (CV/AB) contributes similarly to the soft flux 
of the gas based on Lx~Lk scaling relations within 4" 
for this component, where Lk is the Two Micron All Sky 
Survey (2MASS) if -band luminosity. Hence, including 
this component is essential in the analysis. 

For the LMXB component, we take Tlmxb = I.6I+002 
measured from the combined spectrum of all the re- 
solved point sources within D25 of NGC 3115. For 
the CV/AB component, we fitted the unresolved X-ray 
emission of M32, which is believed to be tru ly gas free 
(jRevnivtsev et al.l 120071 iBoroson et "ail 120111 ). using an 
absorbed thermal + power-law model. The best-fit tem- 
perature is T cv /ab = 0.59^^13 keV and the power-law 
index is r cv /AB = l-93l ; 09 , which are consistent with 
the values measured by iRevnivtsev et all ([2008) . The 
CV/AB normalizations of each annulus are determined 
by the Lx~Lk scaling relation derived from M32. 

We used the X-ray Spectral Fitting PackageQ (XSPEC) 
for spectral analysis. The 2001 and 2010 spectra of 
each annulus in the 0.5-6.0 keV energy range were fitted 
jointly to the three-component absorbed model (PHABS 
in XSPEC)— a thermal model (APEC in XSPEC) for 
the gas, a power-law (POWERLAW in XSPEC) for 
the unresolved LMXBs, and a combination of thermal 
+ power law (APEC + POWERLAW) model for the 
CV/ABs. The absorption is fixed at the Galactic value 
of N H = 4.32 x 10 20 cm" 2 ()Dickev fc Lockmanl fl990h . 
We fit the temperature of the thermal gas component. 
The normalizations of the thermal gas and the LMXB 
components are allowed to vary and untied. Since the 
abundance cannot be constrained from the spectra, we 
fixed it to solar value. Using 0.2 solar only increases 
the best-fit temperature by < 10% and increases the gas 
normalizations by a factor of four. The derived density 

4 http:/ /heasarc. nasa.gov/xanadu/xspec/ 
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Fig. 2. — Upper panel: projected temperature profile of the hot 
gas component. Different colors are with different binnings. Lower 
panel: deprojected density profile. The best-fit power-law model 
within 4" is shown in solid line. 



only increases by a factor of two and all our conclusions 
remain the same. All the other parameters are fixed as 
described above. The spectra were fitted using the c- 
statistic. Varying IYmxb or Tcv/ab by their la uncer- 
tainties or varying the CV/AB normalizations by ±20% 
(that is larger than their la uncertainties) only changes 
the best-fit gas temperature or the gas normalizations by 
<9%. 

3. TEMPERATURE PROFILE 

The projected temperature profile is shown in the up- 
per panel of Figure [21 Although the error bars are quite 
large, it can be seen that the trend is rising toward the 
center as expected in a Bondi flow. We have also con- 
structed a physical model with a temperature profile ex- 
pected in a Bondi flow (i.e., T oc and a density 
profile taken from our best- fit model (p(R) oc see 
Section [4] below) . The simulated projected temperature 
profile is consistent with our measured profile. In par- 
ticular, the simulated projected temperature within the 
2" region is biased low with respect to the deprojected 
(physical) temperature and is consistent with our mea- 
sured low value. Thus, the bias in the central temper- 
ature may be explained by the projection of the cooler 
gas from the outer region. 

For a measured ISM temperature of 0.3 keV in the 
outer most bin, this corresponds to Rb =112-224 pc 
(2'.'4-4'/8) depending on the adopted black hole mass es- 
timate (1-2 x 10 9 Af Q ). Careful selection of the boundary 
radius suggests that the onset in the increase in temper- 
ature begins at 4"-5", rather than at 2". To confirm 
this, we have performed a Markov chain Monte Carlo 
(MCMC) simulation to address the significance of the 
rise in temperature just within a radius of 4". The dis- 
tribution (normalized probability) and the cumulative 
probability of the simulated temperature ratios between 
the 2"~4" and the 4"-10" annuli are plotted in Figured 
The MCMC test suggests that the temperature ratio is 
larger than 1.13 at the 90% confidence level. Assum- 
ing that the temperature beyond Rb is constant and 




Fig. 3. — Normalized probability (solid line) and the cumula- 
tive probability (dotted line) of the simulated temperature ratios 
between the 2"-4" (Ti) and the 4"-i0" (T 2 ) annuli. The vertical 
line cuts the cumulative probability at f0%. 

T oc R^ 1 inside i?_s0, the expected temperature ratio 
at 3" and at Rb (4") is 1.33, which is consistent with 
the peak ratio of 1.4 in the MCMC simulation. 

The harder spectrum within 4" is not the spurious re- 
sult of any hard central point source. We obtained the 
same overall result even when we excised the inner 1". 
The effect of unresolved X-ray binaries has been prop- 
erly accounted for by quantifying the unresolved emis- 
sion in the 2.0-6.0 keV band and extrapolating to lower 
energies. The effect of a large number of softer accret- 
ing white dwarfs, stellar coronae, etc., has also been 
taken into account in our analysis. Such an increase 
in the temperature is exactly what is expected in all 
models in which the black ho le is gravitationally captur- 
ing the ambient gas (see, e.g . , iQuataert fc Naravanll 19991 : 
IQ.uataert fc Gruzinovl [20001 ) . In fact, a similar spike in 
ISM temperature was found at the center of NGC 4649 
and was interpreted as confi rmation of a >10 9 -M© black 
hole (Humphr ey et al.ll2008| ). although in this galaxy Rb 
is unresolved by Chandra. Finally, the lack of detectable 
radio emission from the core of NGC 3115 argues against 
heating of the gas from jets. 

4. SURFACE BRIGHTNESS AND DENSITY PROFILES 

We have calculated the surface brightness profiles for 
the thermal gas, CV/AB, and LMXB components. The 
CV/AB profile was calculated using the same Lx~Lk 
scaling relations mentioned above. The LMXB pro- 
file was calculated by assuming that the emission above 
2 keV is contributed by the LMXBs only, and its contri- 
bution to the soft emission was calculated by the assumed 
spectral model above. The gas profile in the soft band 
was calculated by subtracting the CV/AB and LMXB 
contributions. Figure [4] shows the surface brightness pro- 
files for the three components in the 0.5-1.0 keV energy 
band. The upper limit is chosen to maximize the con- 
tribution from the gas component with respect to the 
other two components. Note that the gas component is 

5 Note that the temperature profile can be flatter in the tran- 
sition region between the ambient IS M and the central accretion 
flow near ~ 0.1-1 R B (Quatacrt 2002). 
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Fig. 4. — Surface brightness profile for the hot gas (black circles), 
CV/AB (red open squares), and LMXB (green triangles) compo- 
nents in the 0.5-1.0 keV band. 

clearly detected within ~ 15" and is rising more steeply 
compared to the stellar (CV/AB) component toward the 
center. 

With the surface brightness profile for the hot gas com- 
ponent, the electron density profile can be deprojected. 
We converted the count rate in each annulus into flux 
using the XSPEC APEC model with a temperature of 
0.45 keV and solar abundance. The temperature was 
taken as the mid value of the best-fit temperature of the 
central 0"-2" bin and the ambient temperature. Since 
the emissivity mainly depends on density, the precise 
value of temperature is not important at our level of ac- 
curacy. We have tested the density profiles derived using 
different temperatures between 0.3-0.6 keV. All the de- 
rived density profiles are consistent within the error bars. 
The flux was then deprojected into e missivity, e, a t each 
radius using a meth od described in | Wong et al.l (2008) 
and also described in lKriss et al.l (|1983| ). This technique 
calculates the emission measures of each spherical shell 
starting from the outer most annulus toward the center, 
and the emission measure of each subsequent shell is cal- 
culated by subtracting the projected emission from the 
outer shells. The electron density can then be converted 
from the emissivity by the relation e = n 2 A, where the 
emissivity function A depends on the temperature and 
the abundance. The emissivity function is taken from 
the same APEC model. 

The deprojected electron density profile is shown in the 
lower panel of Figure [2] The errors were estimated by 
running 10 6 Monte Carlo simulations. Fitting the density 

_r-i qq + 0.231 

profile within 4 to a power law gives p cx R 1 - u -0.21J. 
This is consistent with the power law index of 1.14+0 
in the 4"-25" region. 

5. DISCUSSION 

Accretion models predict that gas within the Bondi ra- 
dius flows toward the SMBH. The increases in the tem- 
perature and density of thermal gas toward the center 
are direct evidence that the gas is being gravitationally 
captured by the SMBH. Resolving the temperature and 
density profiles within the Bondi radius can provide tight 
constraints to accretion models. 



Our 150 ks Chandra observation of NGC 3115 supports 
that there is a rise in temperature inside Rb compared 
to the galaxy's ISM temperature outside Rb, exactly as 
expected from accretion models. We have shown that the 
onset in the increase in temperature begins at » 4"-5" 
rather than at 2". This supports that the mass of the 
SMBH is at the upper end of the 1-2 x 10 9 M W range 
(jKormendv et al.lll996t lEmsellem et al.lll999f ). 

In previous studies where the Bondi accretion flow 
could not be spatially resolved, the overall spectrum 
has to be modeled to constrain accretion model. With 
the spatially resolved spectra we have measured, we can 
draw definite conclusion for the first time on the dy- 
namical properties within the Bondi radius. In these 
regions (R > 100R S , where R s = 2GM B n/c 2 is the 
Schwarzschild radius), the electron temperature is fairly 
well determin ed since the gas is essential ly virial and one 
temperature (jQuataert fc Na ravan 1999). Therefore, ob- 
servations of the bremsstrahlung emission in X-ray give 
direct information on the density of the outer regions 
of the Bondi flow and thereby the accretion rate on the 
outside (M(R > 10Q-R S )). 

Because the uncertainties in our measurements are still 
large, we do not attempt to compare the observational 
results to state-of-art theoretical models. Instead, we 
focus on the simplest c lassical models. In the early an- 
alytic models (ADA Fs; Ilchimardfl977t iRees et alJl'l982t 
iNaravan fc Yilll994[ ) it was argued that the radiative ef- 
ficiency is low because the accretion flow is a hot two- 
temperature plasma in which most of the energy is car- 
ried into the black hole by the ions. If the electrons — 
which produce the radiation we observe — receive very lit- 
tle of the energy in the system, the radiative efficiency 
will be much lower than the standard value of ~10%. In 
ADAF models, the accretion rate is argued to be simi- 
lar to the Bondi rate and the dynamics of the inflowing 
gas is similar to that of spherical Bondi accretion, even 
though the inflow has angular momentum. Contrary to 
ADAF and Bondi predictions, more recent analytic mod- 
els suggest that when the radiative efficiency is low, very 
little of the mass available at large radii actually accretes 
onto the black hole, with most o f it being blown away at 
all radii in an outflow (ADIOS; iBlandford fc Begelmanl 
I1999T ) or co ntinuously circulating in convective ed- 
dies (CDAF; INaravan et alJl2000t iQuataert fc Gruzinovl 
2000). Accretion rates <C Mb are strongly supported by 
global numerical simulations of thick disk accretion (e.g. , 
I lgumenshchev fc Abramowicz 1 119991: iStone et al.l 119991 : 
lHawlev fc Balbusll2002l : llgumenshchev et al . 2003). 

Here, we assume the self-similar structure of the accre- 
tion models extends from th e vicinity of the SMBH ou t to 
the Bondi radius (see, e.g.. lAbramowicz et al.ll2002l ). In 

the classic Bondi/ADAF models, p(R) cx R~ 3 ^ 2 , while in 
CDAFs the redistribution of gas via convection predicts a 
much flatter density profile, p(R) cx R- 1 / 2 . The ADIOS 
models allow the accretion rate to vary with radius and 
the density profile is modified as p(R) cx R~ 3 / 2+p , where 
p w 0-1 characterize the density and accretion rate sup- 
pression. In our 150 ks observation, we have constrained 
p = 0.471q23 from the density profile, which is consistent 
with the value of p « 0.6 measure d from the spatially un - 
resolved spectrum of Sgr A* by iBaganoff et al.l (2003) , 
although they cannot put constraints from their data. 
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Taken at face value, we have barely ruled out the classi- 
cal Bondi/ AD AF model and the classical CDAF model 
at about 2a, provided that the density profile has a single 
power-law index within Rb- 

However, theoretical models also suggest that the re- 
gions we have resolved are still i n transition from the am- 
bient ISM to the accretion flow (|Quataertll2002D and the 
slope in these regions can be different from the asymp- 
totic behavior. The measured density slope within Rb 
is consistent with that beyond. Thus, the ADAF/CDAF 
models might not be ruled out. More theoretical work 
is needed to interpret the data. Nevertheless, the data 
suggest that the transition does not finish until at least 
~ 0.2Rb for either the classical Bondi/ ADAF or the 
CDAF models. The density profile around ~ 0.1-1 Rb 
does provide some direct constraints to theoretical mod- 
els, e.g., thermal conduction can change the power- law 
index of the density by abo ut 0.5 around those regions 
(| Johnson fc Quataert l)2007t) . 

We determined that the Bondi accretion rate of NGC 
3115 is M B = 2.2 x 10~ 2 M© yr" 1 . Assuming a 
10% radiative efficiency, the accretion luminosity at 
Rb (10 44 erg s _1 ) is about six orders of magnitude 
higher than the upper limit of the X-ray luminosity 



of the nucleus. Such a discrepancy can be explained 
if the accretion rate is lower than the classical Bondi 
rate. For example, the accretion rate for the ADIOS 
or CDAF models is suppressed and scales with radius 
as M acc ~ aMs(R/ Rb) p , where a is the dimension- 
less viscosity parameter in the standard thin dis k model 
(|Shakura fc Sunvaevl[T97l iBaganoff et al.1 120031 ). Using 
a typical value of a ~ 0.1 to calculate M acc at Rs, p 
should be close to 1 to explain the discrepancy. It is 
also quite possible that the gas has non-negligible angu- 
lar momentum that should be more important at smaller 
radii. This can reduce the accretion rate near the event 
horizon (e.g., Proga & Begelman 2003). Alternatively, 
the radiation effi ciency can be lower than the canonical 
value of 0.1 (e.g., iHcl 120081 ). 

Currently, the uncertainties are limited by the statis- 
tics of our data. A deep Chandra observation will im- 
prove the statistics to provide useful constraints to ac- 
cretion models. 



We thank Eliot Quataert and the referee for helpful 
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